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5a ; Abstract 

The excess entropy, S e , denned as the difference between the entropies of the liquid and the 
ideal gas under identical density and temperature conditions, is shown to be the critical quantity 
connecting the structural, diffusional and density anomalies in water-like liquids. Based on simu- 
lations of silica and the two-scale ramp liquids, water- like density and diffusional anomalies can be 
seen as consequences of a characteristic non-monotonic density dependence of S e . The relationship 



pair correlation approximation to S e . 



between excess entropy, the order metrics and the structural anomaly can be understood using a 
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The behaviour of water is anomalous when compared to simple liquids for which the struc- 
ture and dynamics is dominated by strong, essentially isotropic, short-range repulsions 
For example, over certain ranges of temperature and pressure, the density of water increases 
with temperature under isobaric conditions (density anomaly) while the self-diffusivity in- 
creases with density under isothermal conditions (diffusional anomaly). Experiments as 
well as simulations suggest that the anomalous thermodynamic and kinetic properties of 
water are due to the fluctuating, three-dimensional, locally tetrahedral hydrogen-bonded 
network. Water-like anomalies are seen in other tetrahedral network-forming liquids, such 
as silica^as well as in model liquids with isotropic core-softened or two-scale pair potentials 

In the case of liquids such as water and silica, a quantitative connection between the 
structure of the tetrahedral network and the macroscopic density or temperature variables 
can be made by introducing order metrics to gauge the type as well as the extent of structural 
order Jf3,Q|. The local tetrahedral order parameter, q te t, associated with an atom i (e.g. Si 
atom in SiC^) is defined as 

ft- = l-|E E (cos^ + 1/3) 2 (1) 
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where ipjk is the angle between the bond vectors ry and where j and k label the four 
nearest neighbour atoms of the same type f|. The translational order parameter, r, measures 
the extent of pair correlations present in the system and is defined as 

r =(!/&) /%(£)- l|de (2) 



where £ = rp 1//3 , r is the pair separation and £ c is a suitably chosen cut-off distance |1S 
Since r increases as the random close-packing limit is approached, it can be regarded as 
measuring the degree of density ordering. At a given temperature, qtet will show a maximum 
and r will show a minimum as a function of density; the loci of these extrema in the 
order define a structurally anomalous region in the density-temperature (pT) plane. Within 
this structurally anomalous region, the tetrahedral and translational order parameters are 
found to be strongly correlated. The region of the density anomaly, where (dp/dT)p > 
0, is bounded by the structurally anomalous region. The diffusionally anomalous region 
((dD/dp)r > 0) closely follows the boundaries of the structurally anomalous region, specially 
at low temperatures. In water, the structurally anomalous region encloses the region of 
anomalous diffusivity while this is reversed in silica. 



The pattern of nested anomalies seen in tetrahedral liquids can be reproduced by a model 
liquid with an isotropic, two-scale ramp (2SRP) pair potential with the crucial difference 
that q te t must be replaced by an icosahedral order parameter, q ico , which is defined for a 
particle i as 
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where Yi m (9,<p) denotes the spherical harmonics of order I = 6 averaged over bonds connect- 
ing particle i with its 12 nearest neighbours. The 2SRP potential has two length scales: 
the hard-core and soft-core diameters, o"o and a% respectively with cro/ai = 0.568. A linear 
repulsive ramp potential, u{r) = (£/i/cri)(<7i — r) acts when the interparticle distance r lies 
between <jq and o\. At low temperatures, the 2SRP liquid shows a maximum in q ico and a 
minimum in r as a function of p, resulting in a well-defined structurally anomalous region. 

Based on simulations of silica and the two-scale ramp(2SRP) liquid, here we show that 
the excess entropy, S e , is the critical quantity connecting the structural, diffusional and 
density anomalies in water-like liquids where S e is the difference between the entropy of the 
liquid and that of the ideal gas under identical density and temperature conditions. 

Molecular dynamics (MD) simulations in the canonical (NVT) ensemble have been per- 
formed for the liquid phase of the BKS model of silica p, Q| for which computational 
details are given in ref-jl^. NVT-MD simulations for the Lennard- Jones (LJ) and NVT 
Monte Carlo simulations for the 2SRP liquid were performed using a 256-particle cubic sim- 
ulation cell. Excess entropy for these systems was estimated to better than 90% accuracy 
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using the pair correlation contribution, S2 16, 
isotropic pair interactions, SyjNkn = —2-np J [g(r) In g(r) — g(r) + 1] r 2 dr where g(r) is the 
pair correlation function |l7|. For a binary liquid mixture of two species A and B interacting 
via isotropic potentials, this expression can be suitably generalised and applied to BKS 
silica to give a total entropy, S = Sid + S2 which is consistently 6.5±0.3% higher than the 
values in ref. \^ for the 3.0 g cm -3 isochore. 

We first consider the connection between S2 and the density anomaly. Figure 1 shows 
the behaviour of the excess entropy, S2, and the configurational energy, (U), as a function 
of density for silica and for the two-scale ramp liquid. The S2 (p) curves of the two systems 
are qualitatively the same even though the dependence of (U) on p is completely different. 
The non-monotonic behaviour of S2 with a well-defined maximum and minimum results in 
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extrema in the total entropy, S = Sid + S e , since Sid is a monotonic function of density along 
an isotherm. At these extrema, (dS/dp)T = V 2 (o./k) where a is the isobaric expansion 
coefficient and k is the isothermal compressibility. Provided the liquid is mechanically 
stable (k > 0), maxima and minima in S must correspond to state points which lie at the 
boundaries of the region of the density anomaly where a equals zero. 

A relationship between S e and a suitably scaled diffusivity, D*, is expected on the basis 



of corresponding states relationships of the form D = Aexp(aS e ) for dense liquids 
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231 ] . We have scaled diffusivity for the binary Si0 2 system using an Enskog-type scaling 



|22| and the resulting D*(p) curves have maxima and minima who locations coincide with 
those of So (p) (cf. Figs. 1(a) and 2(a)). The 2SRP system, for which D(p) data is given in 



refs. 
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reveals the same correlation between 5*2 and D. Figure 2(b) shows the strong 



linear dependence of InD* on S2 with a slope of 1.5 for liquid silica. Consequently, the non- 
monotonic behaviour of S 2 as a function of p at low temperatures should give rise to a region 
of anomalous diffusivity in the pT-plane, even if the underlying diffusional mechanisms in 
simple and tetrahedral liquids are widely different 24], ^ | . 

In the case of supercooled water, the configurational entropy (S con f), defined as the 
number of inherent structure basins accessible to the system, has been shown to mirror the 
variation in diffusivity (InD) with temperature in an analogous manner to the S e variation 
discussed above j^fj. S e is a more convenient quantity to deal with it since it can be readily 
obtained from calorimetric data or simulations and is well approximated by a structural 
measure such as S2. While S e can be defined for any system, S con f can be unambiguously 
defined only for supercooled liquids where there is a clear separation of time-scales between 
interbasin and intrabasin motions. 

From the definitions of 5*2 and r, both of which depend on deviations of g(r) from unity, 
it is to be expected that variations in S2 and r are anticorrelated (cf. Figs. 1(c) and 3(a)). 
Thus, the existence of a maximum in S^p) at high densities implies a minimum in r as a 
function of density and the locus of such minima in the pT-plane is sufficient to determine 
the high-temperature boundary of the structurally anomalous region. It is interesting to 
consider the structural changes that give rise to qualitatively similar behaviour in r(p) (or 
S2 (p)) curves in systems as dissimilar as silica and the 2SRP liquids. In the low p and low T 
regime for silica and other tetrahedral liquids, increasing p tends to destroy local tetrahedral 
order thereby reducing q tet and destructure g(r) thereby reducing r and enhancing S2. Since 



the system is very far from the close-packing limit, density-induced ordering is negligible. 
Only after reaching a critical density when excluded volume effects dominate over local 
tetrahedral order, does further increase in density promote pair correlations and lead to a 
decrease in S2 with increasing p. In the case of the 2SRP liquid, there is a single length 
scale at very low densities, corresponding to the soft wall diameter, <j\ — 1, as shown by 
the g{r) curve at p = 1.0 in Figure 3(b). As density increases the second length scale <7o 
emerges. This leads to a peak at r = 0.568 and the greater structure in g(r) tends to 
decrease S , 2 . Simultaneously, the introduction of additional distances for second and higher- 
order neighbours implies that the g{r) function becomes less structured and plateau-like at 
intermediate distances. The net effect is a maximum in r at p = 1.3, followed by a decrease 
until a minimum value is reached at p—1.9. Further increase in p leads to a more structured 
g(r) as the second short length scale (00) becomes dominant and r begins to increase. 

We now consider the chain of arguments that relates S2 and local orientational order, q. 
£2 and r contain the same structural information; therefore it is sufficient to consider the 
relationship between r and q te t or q ico . The characteristic feature of the water-like liquids 
studied here is that the correlation between r and q breaks down with density because the 
energetically favourable local structures have very different symmetries in the low density 
and high density limits. Systems, such as LJ, characterized by strong, short-range isotropic 
repulsions have local icosahedral order in the liquid phase and hep or fee packing in the 
solid phase. Increasing density enhances particle correlations and local structure leading to 
a strong correlation between r and qi co [12|. In contrast, the stable crystalline forms of the 
2SRP system at low and high densities are fee and rhombohedral respectively suggesting 
that the low density liquid will have local, icosahedral order which will be attenuated with 
increasing density; consequently the correlation between qi CO and r will break down with 
density. In the case of tetrahedral liquids, the local geometry at low and high densities is 
tetrahedral and icosahedral respectively and therefore there is a density-driven breakdown 
of the correlation between q te t and r. Figure 4 illustrates how 5*2 shifts from being negatively 
correlated with qtet/ Qico in the anomalous regime to positively correlated in the normal regime 
along various isotherms as a consequence of the changing relationship between translational 
and orientational order. 

It is useful to apply the above arguments to the one-scale ramp (1SRP) liquid with 
o"o = which displays density and diffusional anomalies but not the structural anomaly 



since qi C o(p) does not have a maximum |lOj ]. The 2SRP and 1SRP liquids show the same 
qualitative behaviour in r and 5*2 with an effective second length scale originating from a 
peak in g(r) close to zero at intermediate and high densities. However, there is no density- 
induced shift in symmetry of local order since the low and high density crystalline structures 
are fee and hep respectively, both of which are compatible with icosahedral local order in 
the liquid. Consequently, 1SRP does not display a well-defined structural anomaly. 

We now consider the relationship between the positional or configurational contribution 
to the entropy (S e ) and internal energy {{U)). Figure 5 shows the striking contrast between 
the S e versus (U) plots for silica, 2SRP and Lennard- Jones (L J) liquids, generated by varying 
density along different isotherms. The LJ system shows a strong positive correlation between 
^2 and (U). In the case of the 2SRP liquid, the positive correlation between S 2 and (U) 
in the anomalous regime shifts to a negative correlation at high densities. The 5*2 versus 
(U) curves for the tetrahedral liquid show distinct low-density and high-density segments. 
In the anomalous regime, the entropy of the high density segment always lies above that 
of the low-density one, while in the normal regime this order is reversed; the two segments 
coincide for the crossover temperature. Since it is straightforward to evaluate S 2 and (U) 
from simulations or experimental data, such plots should provide a convenient means to 
diagnose anomalous behaviour in a range of liquids, including ionic and intermetallic melts 
and complex fluids with ultra-soft repulsions. 

Our results demonstrate that if the excess entropy, S e , of a liquid shows a well-defined 
minimum followed by a maximum as a function of density, then the system will have nested 
diffusional and density anomalies. Since S e is well approximated by the pair correlation 
entropy, S 2 , a connection with the order metrics can be readily formulated to show that if 
the energetically favourable local geometries in the low and high density limits have different 
symmetries, then a structurally anomalous regime can be defined in terms of an orientational 
local order parameter and a translational or pair-correlation order parameter. Such anoma- 
lous liquids will show a characterstic fingerprint when excess entropy is plotted as a function 
of configurational energy. The relationship between the excess entropy and the diffusional 
and density anomalies will be valid even if the pair correlation approximation is substituted 
by more accurate experimental or computational estimates. Provided pair correlations play 
a dominant role in liquid state structure and dynamics, as is likely in most systems of inter- 
est, the close correspondence between S e and the translational order parameter will remain 



and our conclusions regarding the structural anomaly will also be valid. 
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Figure Captions 

1. Dependence of pair correlation entropy, S 2 , and configurational energy, (U), on density, 
p. Results for silica are shown in the left panel with isotherms corresponding from top 
to bottom to T = 6000K, 5500K, 5000K, 4500K and 4000K. The units of p, S 2 and 
(U) are g cm -3 , J K _1 mol" 1 and 10 3 kJ mol -1 . Results for 2SRP liquid are shown in 
the right panel with isotherms corresponding from top to bottom to T = 0.2, 0.109, 
0.082, 0.063, 0.045, 0.036 and 0.027. Reduced units of p, S 2 and (U) are af 3 , k B and 

2. Scaled diffusivity, D*, of liquid BKS silica: (a) Dependence of D*, on p at differ- 
ent T and (b) Correlation plot of D* and S 2 . The scaled diffusivity, D* equals 
(Dsi/xsi) XSi (Do/xo) x ° i where D a , x a and Xa are the self-diffusivity, mole fraction 



and scaling parameter respectively of component a in the mixture [22j. The straight 
line in part(b) corresponds to D — 0.16exp(1.5S , 2); the three lowest density state 
points at 4000K have been excluded from the fit. 

3. (a) Dependence of translational order parameter, r, on p for 2SRP liquids along 
isotherms at T values given in Figure 1. (b) Pair correlation function, g(r), at different 
densities along the T=0.027 isotherm. 

4. Correlation between the excess entropy, S 2 and (a) the tetrahedral order parameter, 
qtet, for silica and (b) the icosahedral order parameter, q ico , for 2SRP liquid. Isotherms 
are shown in order of increasing T from left to right at T values given in Figure 1. 



5. Correlation plot of excess entropy, S 2 and configurational energy, (U) for (a) silica (b) 
2SRP and (c) Lennard- Jones (LJ) along different isotherms. The highest and lowest 
density state points for each isotherm are marked by horizontal and vertical arrows. 
Isotherms for 2SRP liquid labelled as in Figure 1. 
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